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Abstract Isoradial graphs are a natural generalization of regular graphs which give, 
for many models of statistical mechanics, the right framework for studying mod- 
els at criticality. In this survey paper, we first explain how isoradial graphs natu- 
rally arise in two approaches used by physicists: transfer matrices and conformal 
field theory. This leads us to the fact that isoradial graphs provide a natural setting 
for discrete complex analysis, to which we dedicate one section. Then, we give an 
overview of explicit results obtained for different models of statistical mechanics 
defined on such graphs: the critical dimer model when the underlying graph is bi- 
partite, the 2-dimensional critical Ising model, random walk and spanning trees and 
the g'-state Potts model. 



1 Introduction 

Statistical mechanics aims at describing large scale properties of physics systems 
based on models which specify interactions on a microscopic level. In this setting, 
physics systems are modelled by random configurations of graphs embedded in a 
fif-dimensional space. Since vertices typically represent atoms, the goal is to let the 
mesh size tend to zero and rigorously understand the limiting behavior of the system. 

Although real world suggests that we should focus on the case of 3-dimensional 
systems, we restrict ourselves to the case of dimension 2. There are two main rea- 
sons guiding this choice: first, for models we consider, only very few rigorous results 
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exist in dimension 3, and more importantly, 2-dimensional systems exhibit beautiful 
and rich behaviors which are strongly related to this choice of dimension. 

Historically, the most studied graph is certainly I?, followed by the honeycomb 
and triangular lattices. Since solving statistical mechanics models involves dealing 
with combinatorial and geometric features of the underlying graph, it has been most 
convenient to handle the simplest and most regular ones. However in this paper, we 
deal with models defined on a more general class of graphs, called isoradial graphs. 

The motivation behind this generalization is to find the most appropriate setting, 
which exhibits some essential features required to solve the questions addressed, 
thus allowing for proofs revealing the true nature of the problems, and hopefully 
a full understanding of the issues at stake. Before going any further, let us define 
isoradial graphs and the corresponding rhombus graph. 

Definition 1. A graph G = {V,E) is said to be isoradial IIKen02bl . if it has an em- 
bedding in the plane such that every face is inscribed in a circle of radius 1, and such 
that all circumcenters of the faces are in the closure of the faces. 

From now on, when we speak of the graph G, we mean the graph together with a 
particular isoradial embedding in the plane. Examples of isoradial graphs are given 
in Figure [T]below, in particular Z^, the honeycomb and triangular lattices are isora- 
dial. 




Fig. 1 Examples of isoradial graphs: the square lattice (left), the honeycomb lattice (center), and a 
more generic one (right). Every face is inscribed in a circle of radius 1, represented in dashed lines. 



To such a graph is naturally associated the diamond graph, denoted by G*, de- 
fined as follows. Vertices of G* consist in the vertices of G, and the circumcenters 
of the faces of G. The circumcenter of each face is then joined to all vertices which 
are on the boundary of this face, see Figure |2] Since G is isoradial, all faces of G° 
are side-length- 1 rhombi. Moreover, each edge e of G is the diagonal of exactly 
one rhombus Rg of G*; we let 9e be the half-angle of the rhombus at the vertex it 
has in common with e. For later purposes, we label vertices of the rhombus by 
vi , V2, V3, V4, as in Figured (right). 

Isoradial graphs naturally arise in two different approaches to statistical mechan- 
ics. We describe them in the next two sections. 
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Fig. 2 An isoradial graph (left). The white dots are the circumcenters of the faces, which are also 
the vertices of the dual G* . Its diamond graph is represented in the center On the right is the 
rhombus Re and the half-rhombus angle 9(. assigned to an edge e. 



1.1 Transfer matrices, star transformations and Z-invariance 

The first question addressed when solving a model of statistical mechanics is to 
compute the free energy, which is the growth rate of the partition function, count- 
ing the weighted number of configurations. An important technique introduced by 
physicists to solve two-dimensional models is the use of transfer matrices, which 
appears for the first time in the work of Kramers and Wannier IIKW4 1 al IKW4 1 bl . If 
a model of statistical mechanics with local interactions is defined on a torus of size 
n X m, then the partition function Z„j „ can be expressed as the trace of the «-fold 
product of a transfer matrix T, whose rows and columns are indexed by the config- 
urations of the model in a strip mx I, and if ^ and are two such configurations, 
then the matrix element -^z is the contribution to the Boltzmann weight of the 
interactions in "t^ and are in consecutive strips. The free energy of the model / 
can be expressed in terms of the largest eigenvalue of T denoted by Aj ; 

/ = - lim — log(Z,„.„) = - lim — log(Ai). 

If the interaction constants on the graph are not homogeneous, then the trans- 
fer matrices used to pass from one strip to another are different. If we want to be 
able to compute the free energy, as in the homogeneous case, in terms of spectral 
characteristics of the transfer matrices, then these should commute. 




m 



Fig. 3 Transfer matrices for an inhomogeneous model on the m x n torus. 
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Since the interactions are local, the transfer matrices are themselves sums of local 
operators R acting on configurations at neighboring sites, and diagrammatically rep- 
resented by one rhombus of a strip. A sufficient condition to ensure the commutation 
of transfer matrices is to demand that these /^-matrices satisfy some algebraic rela- 
tion, called the Yang-Baxter equations, which can be loosely formulated as follows. 
Consider the two ways to tile a hexagon with the same three rhombi on Figure HI 
The transformation from one to the other is called a star-star transformation. The 
model satisfies the Yang-Baxter equations if the sum of Boltzmann weights over lo- 
cal configurations of the left hand-side equals those of the right hand side. This very 
strong constraint gives a set of equations which needs to be satisfied by Boltzmann 
weights, see IIPAY06l[Bax89l for more details. If a non trivial solution is found, we 



Fig. 4 Star-star transformation (the name star referring to the 3-branches stars). 



say that the model is Z-invariant. 

The condition of Z-invariance yields a natural generalization. Suppose that the 
graph embedded on the torus of size m x n is a periodic tiling by rhombi (where 
the period is independent of m and n). Then by performing a sequence of star-star 
moves, this graph can be transformed into big pieces of tilted copies of IIKen02al . 
If boundary effects can be neglected, then it suffices to solve the model on each of 
the copy of pieces of Z? using transfer matrices. Summarizing, if the model is Z- 
invariant, then it can be solved on any periodic graph consisting of rhombi. 

Consider an isoradial graph G and recall that a unique rhombus of the underlying 
diamond graph can be assigned to every edge of G. Then, looking at edges rather 
than rhombi, the star-star transformation becomes a star-triangle transformation, 
see Figure |5]below. 



Fig. 5 Star-triangle transformation for isoradial graphs. 



The definition of Z-invariance naturally extends to this setting, and if the model is 
Z-invariant the transfer matrix approach can be performed on G*, thus explaining the 
occurrence of isoradial graphs and Z-invariance in the context of transfer matrices. 
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Conformal field theory (CFT), introduced by Belavin, Polyakov and Zamolodchikov 
IIBPZ84I (see also IIDFMS97II ) is a theory which aims at describing models at crit- 
icality, supposed not to depend on specificities of the graph, and to be conformally 
invariant. This very strong statement remained largely inaccessible to the mathe- 
matics community (except for a few models, such as dimers on the square lattice 
BKenOOl ) until the introduction of the SLE process by Schramm in 1999. The SLE 
is conformally invariant and conjectured to be the limiting process for many mod- 
els of statistical mechanics, thus filling a huge gap with CFT. For relations between 
CFT and SLE, see for example IIFW03I IBB04I . Several of these conjectures are 
now solved: loop erased random walk and uniform spanning tree (Lawler, Schramm, 
Werner IILSW04I ). site percolation on the triangular lattice (Smirnov BSmiOll ). Ising 
(Smirnov BSmilOI . Chelkak and Smirnov ICSbl ). . . 

It is thus of key importance to have a setting suitable for doing discrete complex 
analysis, and proving convergence to its continuous counterpart. To this purpose, 
isoradial graphs are a perfectly suited object. Indeed, when considering an edge and 
its dual, they consist in the two diagonals of a rhombus, and are thus orthogonal. 
This allows for a natural discretization of the Cauchy-Riemman equation. Since a 
lot of developments have happened in this direction, we dedicate it the next section. 



2 Discrete complex analysis on isoradial graphs 

As mentioned above, isoradial graphs provide a natural framework for a generaliza- 
tion of the construction in the case of the square lattice. The ideas presented here go 
back to Duffin IIDuf68l . and have been developed later by Mercat IMerOll . Kenyon 
IIKen02bl . Chelkak and Smirnov OCSal . Cimasoni OCiml and others. There are other 
possible discretizations of complex analysis on isoradial graphs with remarkable 
properties. An example using the notion of cross-ratios is given in IIBMS05I . 



2.1 Discrete holomorphic and discrete harmonic functions 

Let us start with the notion of discrete holomorphy. In the continuous setting, the 
Cauchy-Riemann equations satisfied by a holomorphic function imply that its partial 
derivatives along two orthogonal unit vectors differ by a factor ;. In the setting of 
isoradial graphs, one takes advantage of the fact that faces of are rhombi, and 
thus have perpendicular diagonals, to write down finite difference equations which 
are discretizations of Cauchy-Riemann equations. 

Let / be a C-valued function on vertices of G*, let e be an edge of G, and Re be 
the corresponding rhombus with the labeling of its vertices vi, V2, V3, V4, given in 
Figure |2] We say that the function / is discrete holomorphic at the rhombus Rg, if: 
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/(V3)-/(V1) ^ /(V4)-/(V2) 
V3-V1 V4-V2 

More generally we define a discrete operator, d : C^'*^*' — > C^^*^', by: 



df{e) = / (v4 - V2) (v3 - VI ) 
where 



/(V4)-/(V2) /(V3)-/(V1)] 4 



V4-V2 V3-V1 



de.yj ^i{vj^i~Vj+i), (1) 

with indices written mod 4. A function / is then said to be discrete holomorphic on 
G", if df = 0. Similarly, we define an operator d by replacing the coefficients dg ^,. 
by their complex conjugate. 

The operators d and d can be extended to 1 -forms (from C^^^) to C^(^*)): let g 
be a C- valued function on edges of G, then for all v G y(G*), 

e=(v,y')~v ' 

The extension of the operator d is again obtained by replacing coefficients by their 
complex conjugate. 

Note that the normalization of the operators d and d we adopt here differs from 
that chosen by Mercat, or by Chelkak and Smirnov. One can recover their normal- 
ization by multiplying our operators by diagonal operators acting on edges. So these 
variations do correspond to the same notion of discrete holomorphy. 

An notion intimately connected to holomorphy is that of harmonicity. The 
Laplace operator 4 on G can be written as the restriction of dd to functions sup- 
ported on gQ. Its action on a function / defined on vertices of G is given by 

VveG, Af{v)^d{df){v) - ^ tan 0„, (/(«)- /(v)). 

A function / is said to be discrete harmonic if Af = 0. This corresponds to choosing 
for every edge e, a conductance equal to tan Qg. 

From the factorisation A = dd, one sees that, like in the continuous case, the re- 
striction of a discrete holomorphic function to G (resp. G*) is harmonic. Conversely, 
if is a harmonic function on G, there exists a harmonic function H* on G* such 
that H + iH* is discrete holomorphic. 

From this choice of Laplacian, one can define the discrete analogues of classical 
quantities in continuous potential theory (Green function, Poisson kernel, harmonic 
measure,. . . ). A surprising amount of estimates have a discrete version, and can be 
found in llCSal . 

Discrete complex analysis can serve as an approximation of the continuous the- 
ory. Mercat IMerOlll proved that the pointwise limit of a sequence of discrete holo- 



One can define in a similar way a Laplacian on G* by restricting the same operator dd ioG* 
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morphic functions on isoradial graphs in a domain with a mesh going to zero is 
holomorphic. See QCiml for a generalisation of this result to isoradial discretization 
of compact manifolds. These local convergence results have been supplemented by 
global convergence theorems OCSbl . implying in particular the uniform convergence 
for the discrete potential theory objects to their continuous counterparts. 



2.2 Discrete exponential functions 

There is a class of functions playing a special role in this theory. These functions, 
called discrete exponential functions, are defined recursively on G" . For any given 
vertex vq of G" and any A £ C, the function Ex^^,^{-\X) is defined as follows: its 
value at vq is 1, and if vi and V2 are neighbors in G*, let e'" ~ V2 — vi, then: 

Exp,„ (V2; A ) = Exp,„ (vi ; A ) • 

The name comes from the fact that these functions satisfy the following identity: 
for any pair of neighbouring vertices vi and V2 of G", 

V f V f 7^ Exp (v2;A)+Exp (vi;A) 

Exp„g (v2; A ) - Exp,,^ (vi ; A ) = A (v2 - vi ) , 

which is a discrete version of the differential equation 

dexp(Ax) = A exp(Ax)dJc. 

satisfied by the usual exponential function. It is straightforward to check from the 
definition that Exp^,^{-;X) is discrete holomorphic. In fact, Bobenko, Mercat and 
Suris IIBMS05II show that any discrete holomorphic function can be written as a 
(generalized) hnear combination of discrete exponential functions, at least under 
the quasicrystallic assumption (a finite number of slopes ±e'^^ , ±e'^'' for the 
rhombi chains): if / is a discrete holomorphic function on G*, then there exists a 
function g{X) defined on a neighborhood in C of {±e'^i , . . . , ±e'^'' } such that 

Vv e G^ /(v) = j( Exp,Jv;A)g(A)dA, 

where F is a collection of disjoint small positively oriented contours around the 
possible poles ±e'^J of the integrand. 



8 Cedric Boutillier and Beatrice de Tiliere 

2.3 Geometric integrability of discrete Cauchy-Riemann equations 



An important feature of the discrete Cauchy-Riemann is the property of 3D consis- 
tency IIBMS05I . Let / be a discrete holomorphic function on G*. Given its value on 
three vertices of a rhombus, the other can be computed using the discrete Cauchy- 
Riemann equation. 

If a star-star transformation is reaUsed on the graph (see Figure |4|i, can we 
determine the value of the function / at the new vertex in the center, so that / is 
again discrete holomorphic? A priori, there are three different ways to compute a 
possible value: one for each rhombus. It turns out that the three values obtained are 
equal. 

Under the quasicrystallic assumption, G" can be seen as a monotone surface in 
l/, projected back to a properly chose plane. The star-star move is a local displace- 
ment of this surface along the faces of a cube. If given two monotonic surfaces I.\ 
and L2 in Z'' that can be deduced one from the other by a sequence of star-star 
moves, then there is a canonical way to extend a discrete holomorphic function / 
defined on Li to L2 by "pushing" its values along the deformation making use of 
the Cauchy-Riemann equations. In this sense, this 3D consistency can be consid- 
ered as a notion of integrability IIBS08II . This property is very much related to the 
Z-invariance in statistical mechanics. 



2.4 Generalization of the operator d 

The double graph G^ of G is a planar bipartite graph constructed as follows, see 
also Figure|6]below. Vertices of G^ are decomposed into two classes: black vertices, 
which are vertices of G and of G*, and white vertices, which are the centers of the 
rhombi of G*, seen either as edges of G or as those of G*. A black vertex b and 
a white vertex w are connected by an edge in G^ if the vertex corresponding to b 
is incident to the edge corresponding to w in either G or G*. That is, edges of G^ 
correspond to half diagonals of rhombi of G*. It turns out that G^ has an isoradial 
embedding (with rhombi of side-length 5) obtained by splitting each rhombus of G* 
into four identical smaller rhombi. 

The operator d : C^'^^"^ -> C^('^) of Section O can be interpreted now as an 
operator d : C^C^") C^ic")^ where B(G^) (resp. W{G°)) denotes the set of black 
(resp. white) vertices of G^. Let / be a C-valued function defined on black vertices 
of G^, then for every white vertex of G^: 

Let w,x,b,y be the vertices of a rhombus enumerated in cclw order, so that w G 
W(G'^) and beBiG°). Then by Formula O, 
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Fig. 6 The double graph of the graph on the right of Figure[T] in thick full lines is the original 
graph G, in full lines is the dual graph G* , in dotted lines is the underlying G* graph. 



This definition can be extended to any bipartite isoradial graph, see Section lTTI 

For a detailed study of this operator on isoradial graphs on compact surfaces and 
its connections with discretization of geometric structures (especially spin struc- 
tures), see llCiml . 



3 Dimer model 

The dimer model represents the adsorption of diatomic molecules on the surface of 
a crystal. The surface of the crystal is modeled by a planar graph G ~ iy^E), which 
we assume to be finite for the moment. A dimer configuration of G is a perfect 
matching of G, that is a subset of edges M such that every vertex of G is incident 
to a unique edge of M. Denote by ^(G) the set of dimer configurations of G. An 
example of dimer configuration is given in Figure |8] Assume that a positive weight 
function V is assigned to edges of G, that is every edge e of G has weight Vj. . Then, 
when the graph G is finite, the probability of occurrence of a dimer configuration M 
is given by the dimer Boltzmann measure: 



where Zdimer = Lmg^(G) rieGM is the normalizing constant known as the parti- 
tion function. 

The dimer model has the attractive feature of being exactly solvable, meaning 
that there is an explicit expression for the partition function. This fundamental result 
is due to Kasteleyn ||Kas61| and independently to Temperley and Fisher IITF61I . It 
relies on the Kasteleyn matrix, which we now define. A Kasteleyn orientation of 
the graph G is an orientation of the edges such that every face is clockwise odd 
meaning that, when traveling clockwise around edges of a face, an odd number 
of them are co-oriented. The corresponding Kasteleyn matrix K is the associated 



dn.h = i{x-y). 



(2) 



^'dimer(M) 
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signed, weighted, adjacency matrix of the graph: hnes and columns of K are indexed 
by vertices of G, and for any two vertices x,y of G, the coefficient K^ y is: 

if X and the edge xy is oriented from x to y 
Kx,y = '{ —Vxy if X ^y and the edge xy is oriented from y to x 
otherwise. 



Theorem 1 ( ||Kas61| ITF61I ). The dimer partition function is, 

Zdimer-|Pf(^)|- Vdet(/:). 

Remark 1. The Pfaffian, denoted by Pf, of a skew-symmetric matrix is a polynomial 
in the entries of the matrix, which is a square root of the determinant. 

Proof (sketch). Refer to IIKas67l for details. When writing out the Pfaffian of an 
adjacency matrix as a sum over permutations, there is exactly one term per dimer 
configuration. The issue is that each term comes with a signature. The purpose of 
the Kasteleyn orientation is to compensate the signature with signs of coefficients. 

As a consequence of Theorem [T] Kenyon derives an explicit expression for the 
dimer Boltzmann measure. 

Theorem 2 ( 0Ken97l ). Let E ~ {ey ~ x\y\ , • • • ,eA- = x^yj^} be a subset of edges of 
the graph G. Then, the probability that these edges are in a dimer configuration 
chosen with respect to the dimer Boltzmann measure is; 



^'dimei({ei,--- ,ek}) 



n^:.,,,, )pf((^: 



where [K ' )e is the submatrix of K ^ whose lines and columns are indexed by 
vertices ofE. 

Proof. Since the proof is very short, we repeat it here. The probability of dimer 
configurations containing edges of E is given by the weighted number of these con- 
figurations divided by the weighted number of all configurations. By expanding the 
Pfaffian along hnes and columns, this yields: 



Using Jacobi's formula for Pfaffians llWOOl , | Pf (/Tfic ) | = | ?f{K) Pf ( {K-'^)e)\ yields 
the result. 

In Section im we state Kenyon's results for the d operator on infinite, bipartite, 
isoradial graphs. Then, in Section [32] we relate them to the corresponding dimer 
model. 
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3.1 Dime operator and its inverse 

As mentioned in Section l24l the d operator, also referred to as the Dime operator 
in IIKen02bl . can be extended to infinite bipartite isoradial graphs. Let G be such a 
graph, the set of vertices can be divided into two subsets B U W, where vertices in 
B (the black ones) are only incident to vertices in W (the white ones). Then d maps 
to C^: let / be a C-valued function on black vertices, then 

(a/)(w)=£a„,,,/(fe), 

where, if w,x,b,y, are the vertices of the rhombus R„i, in cclw order, d„b — K^^y) 
is given by Formula (|2]i. 

One of the main results of the paper IIKen02bl is an explicit expression for the co- 
efficients of the inverse 5, as a contour integral of an integrand which only depends 
on a path joining the two vertices. In order to state the result, let us introduce the 
following functions defined on vertices of G. Fix a reference white vertex wq and let 
A be a complex parameter. Let v be a vertex of G, and wq = vo, , • • • , Vyt = v be a 
path in G" from wq to v. Each edge VjVj+\ has exactly one edge of G (the other is an 
edge of G*). Direct it away from this vertex if it is white, and towards it if it is black, 
and let e"^' be the corresponding vector. The function /,, is defined inductively along 
the path, starting from /v„ = 1 , and 



fy. {X){X — e'"') if the edge vjvj+i leads away from a black or towards 
a white vertex, 



Otherwise. 



It is easy to see that the function is well defined, i.e. independent of the choice of 
path from wq to v. An important point is that, by using a branched cover of the 
plane over wq, the angles aj are defined in M and not [0,27r[. Then, Kenyon has the 
following theorem. 

Theorem 3. Coeffieients of the inverse Dirac operator can explicitly be expressed 
as: 

Cd-\.-o = ^ /}.(A)log(A)dA, (3) 

where C is a closed contour oriented counterclockwise, containing all poles of the 
integrand and with the origin in its exterior. 

Remark 2. The remarkable features of this theorem are the following. Explicit com- 
putations become tractable, since they only involve residues of rational functions. 
Moreover, the Formula (O has the surprising feature of being "local", meaning that 
if the graph is changed away from the vertices wo, b and a path joining them, the 
value of the inverse Dirac operator stays the same. 
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3.2 Dime operator and dimer model 



The Dirac operator d can be represented by an infinite matrix K, whose lines (resp. 
columns) are indexed by white (resp. black) vertices of G, and the coefficient K^^. l, 
is given by: 

K„.b = Kb = iix-y) = 2sin(0„„fc)e'«"", 

where Q„b is the half-angle of the rhombus and e'""* is the unit-vector in the 
direction from w to h. The matrix K resembles a Kasteleyn matrix of the graph G, 
but differs in two instances. 

1. Rows and columns are indexed by only "half" the vertices. To overcome this, 
one can define a weighted adjacency matrix K whose lines and columns are 
indexed by all vertices of the graph: ^ = ( _\t q ) . Then, we would have 

Pf (^) = ± Vdet(£) = ± I A&V{K) I . 

2. The weights of the edges have an extra complex factor of modulus one, instead of 
a sign given by a Kasteleyn orientation. In IIKen02bl . Kenyon shows that the ma- 
trix is Kasteleyn flat, meaning that for each face of G whose vertices are labeled 
by Ml , vi , • • • , M„,, v,„, the coefficients of the matrix K satisfy: 

arg[/:„,,, • ••/:„„„,„,] = arg[(-i • ••/:„„,.„,]. 



Using a result of Kuperberg | Kup98| , this implies that when the graph is finite 



and planar, K behaves as a usual Kasteleyn matrix, thus yielding an explicit for- 
mula for the partition function. 

One of the main conjectures of the paper IIKen02bl is to show that the inverse d 
operator, which we write K^^ using the matrix notation, is related to the dimer 
model on the infinite, bipartite, isoradial, graph G. Its relation to the Gibbs measure 
is given in the following. 

Theorem 4 ( [idT07cl ). There is a unique Gibbs measure -Pdimer dimer config- 
urations of G, such that the probability of occurrence of a subset of edges E{ei = 
w\b\^ - ■ ■ ,6/1 = WkbkS in a dimer configuration ofG chosen with respect to the Gibbs 
measure Pdimer '■S-' 

^dimer({ei,---,ea)= (n^HvA I det [(/T-i )/,, ], (4) 
\i=\ ) 1^'-'^*^ 

where {K-^)i,..ys,j = isgivenby^. 

Proof {sketch). Fix the edge set E, and cut out a finite piece of the rhombus graph 
G" containing E and paths joining vertices of E. Use a result of MdT07cl by which 
any finite piece of the graph G* can be filled with rhombi in order to become part 
of a periodic rhombus tiling of the plane. Define the probability of the edge set E 
as the weak limit of the Boltzmann measures on the natural toroidal exhaustion of 
the periodic graph. Use the uniqueness of the inverse Dirac operator and its locality 
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property to deduce that this expression coincides with (|4|. Use Kolmogorov's ex- 
tension theorem to show existence of a unique measure on G which specifies as (|4|i 
on cylinder events. 



3.3 Other results 

Free energy. Assume that the infinite bipartite isoradial graph is also periodic, and 
let Gi = G/Z^. In |Ken02b|, Kenyon proves an explicit expression for the "loga- 
rithm of the normalized log of the Dirac operator", which has the surprising feature 
of only depending on rhombus angles of Gi : 



Kenyon conjectures it to be related to the free energy of the corresponding dimer 
model, result that we prove this conjecture in ||dT07bl . 

Gaussian free field. There is a natural way of assigning a height function to dimer 
configurations of G. In OdTOTal . we show that when dimer configurations are chosen 
with respect to the Gibbs measure of Theorem|4l the fluctuations of the height func- 
tion are described by a Gaussian Free Field. Note that this proof holds for all dimer 
models defined on bipartite isoradial graphs, in particular for Z? and the honeycomb 
lattice. 



4 Ising Model 

The Ising model, first considered by Lenz IILen20l . has been introduced in the 
physics literature as a model for ferromagnetism. The vertices of a graph G~{V,E) 
represent atoms in a crystal with magnetic moment reduced to one component that 
can take two values ±1. A spin configuration is thus a function a on V, with values 
in {+1,-1}. Magnets with opposite moments tend to repel each other, which has a 
cost in terms of energy. From these considerations, the energy of a configuration a 
on a finite graph is defined by: 



where the positive numbers (Je) are called interaction constants. The probability 
of a occurrence of a spin-configuration is then defined using the Ising Boltzmann 
measure: 




Asing(c^) 
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The normalising factor, Zjsing — Y,a sxp{—S'{a)) is the Ising partition function on G. 

If the graph is drawn on an orientable surface with no boundarjQ, then the parti- 
tion function can be written as a combinatorial sum over contours. There are in fact 
two well-known such expansions: the so-called high-temperature expansion (con- 
tours on G) and the low-temperature expansion (contours on G*). The contours in 
the low-temperature expansion have a nice interpretation in terms of the dual closed 
curves separating zones of different spins. 

The most studied case is the Ising model on (pieces of) regular lattices where the 
interaction constants are taken to be constant equal to j3, the inverse temperature. 
Kramers and Wannier IIKW41all discovered a duality for the Ising model on the 
square lattice: the measure on contours on G* coming from the low-temperature 
expansion is equal to that obtained by considering the high temperature expansion of 
the Ising model on G* at another temperature j3 * satisfying sinh(2j3)sinh(2j3*) = 1. 

This showed that the self-dual temperature for which )3 = j3 * = log y\+\/2 is the 
critical point if this one is unique. Later, Lebowitz and Martin-Lof IILML72I proved 
that this is indeed the case. 

Note that the case of I? is singular since the graph is isomorphic to its dual. 
This duality can be generalized to any graph and its dual. For example the low- 
temperature contour measure for the Ising model on the honeycomb lattice with 
inverse temperature /3 is the same as the high-temperature measure for the Ising 
model on its dual which is the triangular lattice. But using star-triangle transforma- 
tions, one can go back to the Ising model on the honeycomb lattice. This generalizes 
the duahty we had for the square lattice. 

On isoradial graphs, instead of taking the same interaction constants for all edges, 
it is natural to take Je to be a function of the angle Qe- In order to ensure some 
integrability at the discrete level, we can impose that the Boltzmann weights satisfy 
the star-triangle relations, yielding a one-parameter family of Z-invariant interaction 
constants. There is a generalized duality relation: the high temperature expansion on 
G corresponds to the low temperature expansion on G* for a dual value of parameter. 
We qualify as critical the interaction constants corresponding to the self-dual value 
of the parameter, given by the following formula. 



and refer to the corresponding Ising model as critical Z-invariant. See IIBax89L 
chapter 7.13 for the parametrization of the star-triangle relation. They coincide with 
the critical inverse temperatures of the homogeneous triangular, square, and hon- 
eycomb lattices for Oe identically equal to ^, j, | respectively. The Ising model 
on two-dimensional graphs is in correspondence with other well-known models of 
statistical mechanics: the dimer model and the ^-random cluster model with q = 2. 
We differ the discussion of the correspondence with dimers to Subsection l4.2l 



^ The correspondence can be extended to suifaces with boundary by including in addition to the 
closed contours a certain number of paths connected to boundary. 
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There are essentially two ways of taking the large graph limit of the Ising model. 
Either you let the mesh tend to zero at the same time as the number of vertices of 
the graph goes to infinity, in order to get continuous objects in a bounded region of 
the plane (macroscopic level); or you let the graph tend to infinity, keeping the mesh 
size fixed, to get a model of statistical mechanics on an infinite graph. 

The first approach has been adopted by Chelkak and Smirnov OCSbll to prove 
conformal invariance of the Ising model on bounded domains of the plane with 
Dobrushin boundary conditions. The second approach was used by the authors in 
BBdTlObl IBdTlOal to construct a Gibbs measure for the critical Ising model on 
infinite isoradial graphs. 

Before explaining these two approaches, let us mention some works by physicists 
on the Z-invariant model at and off criticaUty by Au-Yang and Perk ||AYP87|[AP04| 
IAP07I and by Martinez ||Mar97l |Mar98l . 



4.1 Conformal invariance 

Let i2 be a bounded domain of M^, and fix two points a and b on the boundary. For 
every £, let G\ be an approximation of 12 by a rhombus-graplH with rhombi of side 
length £. Let also and be approximations of a and b, located at the center of 
half-rhombi on the boundary. 

Consider the Ising model on the isoradial graph corresponding to black ver- 
tices of Gg, with Dobrushin boundary conditions, i.e fixed spins +1 on one arc from 

to b^, and —1 on the other Introduce, following Chelkak and Smirnov ICSbL a 
twisted version of the partition function, considered as a function on edges of G^, 
or equivalently on inner rhombi of G^. Let z denote a generic rhombus of G%, 



ZGe,«E,z =^£(z) = ( sin^ j ^ I Y\ tany 

/ closed contours \e: piece of contour 
+ curve y.ae-^z 



-4wind(y) 



where wind (7) is the winding of the curve y from a to z- If we remove the prefactor 
and the contribution of the winding, this would be for z = be the low temperature ex- 
pansion of the partition function of the Ising model on Qe with Dobrushin boundary 
conditions between Oe and bg. 

Now define for all inner rhombi z of Gg, 

up to a multiplicative factor, introduced for technical reasons. 



^ On the boundary, we put only half-rhombi such that only "black" vertices are exposed on the 
boundary. 
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Chelkak and Smirnov prove that the function Fg and is discrete holomorphic by 
comparing configurations differing by local arrangements near z, and that it solves 
some discrete Riemann-Hilbert boundary value problem. They then deduce using 
approximation results OCSal that, as e 0, the function converges to the function 
solving the analogue continuous Riemann-Hilbert boundary value problem, which 
is conformally covariant. 

Moreover, they prove that this observable satisfies a martingale property with 
respect to the growth of the curve from to bg: given the first n steps of the curve 
(7o = a£,Yu---,Yn), the expected value of the observable ^GE\[yo,%+i],rn+i.fcE W 
equal toFGA[ro,%],r„.''.(2)- 

This martingale property together with the convergence to a conformally covari- 
ant object is sufficient to imply the convergence of the interface between a and b to 
a chordal SLE, with parameter k" = 3. 

In the same paper, Chelkak and Smirnov construct another discrete holomorphic 
observable using not the loops separating clusters of spins, but those separating the 
clusters in the corresponding random cluster modeQ with q^2. This observable is 
a direct generalization of the one introduced by Smirnov in llSmilOII for the square 
lattice. Again, this observable is discrete holomorphic and satisfies a martingale 
property. In this case, the scaHng limit of the interface is a chordal SLE(16/3). 



4.2 The two-dimensional Ising model as a dimer model 

It turns out that a general Ising model on a graph drawn on a surface without bound- 
ary can be reformulated, through its contour expansion, as a dimer model on a dec- 
orated graph. This correspondence is due to Fisher IIFis66ll . Since a lot of exact 
computations can be carried out on the dimer model, this correspondence has been 
proven to be a useful tool to study the Ising model IIMW73II ( IIKas67ll . . . . ) 

The idea is to construct a decorated version of the graph G with a measure 
preserving mapping between polygonal contours of G and dimer configurations of 
. This is an example of holographic reduction II ValOSI . We now present the ver- 
sion used by the authors in IBdTIObl IBdTlOal . which has the advantage of using 
decoration with cyclic symmetry. For other versions of Fisher's correspondence, see 
llFis66l . ICimlOL 

The decorated graph, on which the dimer configurations live, is constructed from 
G as follows. Every vertex of degree ^ of G is replaced by a decoration consisting 
of 3A: vertices: a triangle is attached to every edge incident to this vertex, and these 
triangles are linked by edges in a circular way, see Figure|2]below. This new graph, 
denoted by , is also embedded on the surface without boundary and has vertices 
of degree 3. It is referred to as the Fisher graph of G. 

Here comes the correspondence: to any contour configuration C coming from 
the high-temperature expansion of the Ising model on G, we associate 2l^'*^^l dimer 



* also known as the Fortuin-Kasteleyn percolation. 
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Fig. 7 Left: a vertex of G with its incoming edges. Right: coiTesponding decoration in ; 



configurations on edges present (resp. absent) in C are absent (resp. present) in 
the corresponding dimer configuration of ^. Once the state of these edges is fixed, 
there is, for every decorated vertex, exactly two ways to complete the configuration 
into a dimer configuration. Figure [8]below gives an example in the case where G is 
the square lattice T/? . 



1 2 



3 4 




Fig. 8 Polygonal contour of T? , and corresponding dimer configurations of the Fisher graph. 



In order to have a measure preserving mapping with dimers, the dimer weight 
function v is given for the special case of the Z-invariant critical Ising model by, 

Ve = cot f if e comes from a rhombus of G with half-angle 9, 
Ve = I if e belongs to a decoration. 

A Kasteleyn matrix K is constructed on One of the main result of IBdTlOal 
is that its inverse has a 'local' expression in the spirit of that obtained by Kenyon in 
IIKen02bl for the inverse Dirac operator, see also Theorem[3] 

Theorem 5 ( HBdTlOal ). Ifx (resp. y) is a vertex ofi^ and belongs to the decora- 
tion corresponding to the vertex x (resp. y) of G, then {K^^)x,y has the following 
expression: 

(K-'ky = ^£ g..(Ak,(-A)Exp,(y;A)log(A)dA, 



where g^ (resp. gy) is a simple rational fraction of A, that depends only on the ge- 
ometry of the decorated graph in an immediate neighborhood of x (resp. ofy), and 
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the contour of integration C^.y is a simple closed curve containing all poles of the 
integrand, and avoiding the half-line starting from x in the direction ofy. 

Then, as for the dimer model, we obtain an expHcit expression for a Gibbs mea- 
sure, and recover the expHcit formula for the partition function, obtained by Baxter 
IIBax86l . at the critical point, see LBdTlOaJ for precise statements. 



5 Other models 

We briefly discuss now some aspects of other models of statistical mechanics; first 
some models related to the Laplacian, then the q'-Potts model and its random cluster 
representation, and finally the 6-vertex and 8-vertex models. 



5.1 Random walk and the Green function 

We already mentioned the Laplace operator on functions of vertices, corresponding 
to conductances on edges given by c(e) ~ c{9e) = tanSg. For this particular choice 
of conductances, the associated random walk is a martingale and the covariance 
matrix associated to a step is scalar. Up to a time reparametrization, its scaling limit 
is a standard two-dimensional Brownian motion [CSa Bef08|. Kenyon | K en02b| 
proves a local formula for the inverse of the Laplacian, the Green function 4 ^ ' , on 
an infinite isoradial graph: 



where C is a positively oriented contour containing in its interior all the poles of the 
discrete exponential function, and the cut of the determination of the logarithm. He 
also proves precise asymptotics of G, improved in IBOSI . This Green function gives 
information on the random walk on G with conductances tan 0, but can also be used 
to gather properties of other models from statistical mechanics with an interpretation 
in terms of electric network, such as random spanning trees. 

A spanning tree of G is an acyclic connected subgraph of G whose vertex set 
contains all vertices of G. On a finite graph, if the weight of a spanning tree equals 
the product of the conductances of its edges, the partition function is given, via 
Kirchhoff's matrix-tree theorem IIKir47l . by the determinant of any principal minor 
of the Laplacian. One can construct a measure on spanning trees of G as limit of 
measures on finite graphs where a spanning tree would have a probability propor- 
tional to its weight. The statistics of edges present in the random spanning tree are 
given by a determinantal process on edges; the probability that edges {e\ — v\W\, 
ei^^ VjfWi^} are present is given by 
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where H is the impedance transfer matrix IIBP93II defined by 

H{ei ,ej)=A-^(x'i,Vj)-A-\vi,Wj)-A-\wi,Vj)+A-^ ( w,- , w,- ) • 



5.2 q-Potts models and the random cluster model 

A natural generalization of the Ising model is the Potts model with q states (or 
colors). Spins at neighboring vertices interact only if they have different colors. 
The energy of a spin configuration o is ^(<7) = — Le=i,i G£-^eIo-„^cT„- For q = 2,we 
recover the Ising model. This model can be mapped to the (inhomogeneous) random 
cluster model |Gri06l as follows: 

• spins of different colors are in different percolation clusters, 

• neighboring spins sharing an edge e of the same color are connected with proba- 
bility, pe = 1 -e^^-''. 

Once again, we want the interaction constants to be functions of the half-angle of 
the corresponding rhombus. It is possible to take them to satisfy the Yang-Baxter 
equation. Moreover, if we impose a generalized self-duality, then there is a unique 
choice for the function p ||Ken02a|| : 



Although a formula for the free energy of that model is known ||Bax891 , little is 
known except for the case q = 2 corresponding to the Ising model. 



5.3 6-vertex and 8-vertex models 

To conclude, let us briefly mention that the 6- and 8-vertex models, introduced orig- 
inally on the square lattice, have natural generalization on isoradial graphs, or more 
precisely on the dual of its rhombic graph G**, which is a 4-regular graph. A con- 
figuration of the 6-vertex model (resp. of the 8-vertex model) is an orientation of the 
edges of G"* such that the number of ingoing edge (and thus of outgoing edges) at 
each vertex is equal to 2 (resp. is even). These models can be solved through Bethe 
Ansatz on the square lattice (in the sense that their free energy can be computed 
IILie67llBax7II ). and possess a very rich algebraic structure. There are some studies 
on their Z-invariant generalizations IIBax87llBax86l . but many questions still need 
to solved. 




with r ~ cos {^fqjT) 
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